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^ ■ Abstract 

Ch I We present a theoretical study of the phase diagram and the structure of 

I ' a fluid adsorbed in high-porosity aerogels by means of an integral-equation 

approach combined with the replica formalism. To simulate a realistic gel 
environment, we use an aerogel structure factor obtained from an off-lattice 

C^ . diffusion- limited cluster-cluster aggregation process. The predictions of the 

theory are in qualitative agreement with the experimental results, showing a 
F^ I substantial narrowing of the gas-liquid coexistence curve (compared to that 

^H ' of the bulk fluid), associated with weak changes in the critical density and 

temperature. The influence of the aerogel structure (nontrivial short-range 
correlations due to connectedness, long-range fractal behavior of the silica 
strands) is shown to be important at low fluid densities. 
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I. INTRODUCTION 

The thermodynamics of fluids imbibed in disordered porous materials is a subject of 
great interest, both from an experimental and from a theoretical perspective |l|. Of special 
importance is the question of fluid-fluid phase transitions which are found experimentally 
to be drastically altered compared to the bulk situation. One of these transitions is the 
so-called capillary condensation observed during the adsorption of a single- component gas 

^ . at low temperature. It is marked by a steep increase in the adsorbed mass that usually 

Q i occurs at some pressure below the saturated vapor pressure of the bulk gas. 

From an experimental point of view, one has to distinguish two types of systems. A flrst 
class of adsorbents consists of mesoporous materials of low porosity (between 30 and 60%), 

^ ■ like porous glasses and silica gels. Inside such a matrix, capillary condensation is usually 

accompanied by a pronounced hysteresis loop that shrinks and eventually disappears at some 
temperature Th less than Tc, the bulk critical temperature. This (reproducible) hysteresis 
loop prevents the direct observation of a thermodynamic, equilibrium phase transition, and 
the mere existence of a true phase transition underlying capillary condensation in these 
systems remains a controversial issue. The situation appears clearer for a second class of 
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adsorbents, namely very dilute silica networks like aerogels with porosities as high as 95- 
98% @-§]. There is indeed experimental evidence that fluids in such media undergo true 
phase transitions that terminate at a genuine critical point. Despite the very low density 
of the adsorbing solids, these transitions display features strikingly different from those of 
bulk fluids. Thus, the phase diagrams of ^He and N2 in aerogels show a marked narrowing 
of the liquid-vapor coexistence curve and a displacement of the critical point to a lower 
temperature and a higher density compared to that in the bulk ||^. A similar behavior has 
been observed for a mixture of isobutyric acid and water in a very dilute silica gel ||^ . 

Recently, the study of a disordered lattice-gas model by means of mean-field density 
functional theory has provided new insight into the physical origin of the phenomenology 
associated with capillary condensation |^. First, it has been found that true equilibrium 
phase transitions may occur only when the perturbation induced by the solid is sufficiently 
small. It is thus likely that true phase transitions may occur in aerogels, but not in dense 
matrices of low porosity. Secondly, it has been shown that hysteresis can occur with or 
without an underlying phase transition and that, if a transition exists, the disappearance 
of the hysteresis loop is not associated to capillary criticality. This results from the fact 
that there is a range of temperature and chemical potential for which the system has many 
metastable states ; it relies on the assumption, expected to be valid for low-porosity solids, 
but not necessarily for very dilute gels, that the fluid is not able to fully equilibrate after a 
change in chemical potential (or pressure). 

In the last few years, the concept of "quenched-annealed" binary mixture has been widely 
used to study the structure and the phase behavior of fluids adsorbed in disordered porous 
matrices. In this approach, pioneered by Madden and Glandt and further developed by 
the application of the replica method [0,|[, the fluid molecules equilibrate in a matrix of 
particles frozen in a disordered conflguration that is sampled from a given probability distri- 
bution. On this basis, computer simulations and integral-equation theories using extensions 
of standard liquid-state approximations have been proposed 0. But when the physics of 
adsorption is dominated by the existence of many metastable states (as shown in Ref. 0), 
it is expected that these methods, that have proved very efficient in dealing with bulk flu- 
ids and fluids in simple external potentials, could run into difficulties and overestimate the 
possibility of true phase transitions. 

The case of very dilute silica gels is the most likely exception in which these approaches 
could remain fully operational. Accordingly, we report in the present paper an investigation 
of the phase behavior of a simple fluid adsorbed in an aerogel within the framework of the 
integral-equation theory of quenched-annealed binary mixtures. This study is based on the 
optimized cluster theory (OCT) approach |TD[ whose application to quenched-annealed mix- 



tures was first considered in Ref. |jTT[. It can be considered as an extension of previous work 
toward a description of the actual experimental systems, since, at variance with previous 
studies of quenched-annealed mixtures, a more realistic model of the solid is used. Indeed, 
aerogels display very specific structural properties that may play an important role for fiuid 
adsorption. In the past, this has motivated the development of various lattice models taking 
into account all or part of the specific features of the aerogel- induced disorder ; the phase 
transitions of these models have been studied by various methods, both theoretically and by 



computer simulation [p!^-[T6|. The present work relies on the off-lattice version [0 of the 



model used in Refs. ||T^,|T^ , where the aerogel was described as a connected cluster of sites 



generated by a diffusion limited cluster-cluster aggregation (DLCA) process |T8|JT9|] . 

In the next section, the theoretical framework used in this paper is summarized. In 
section |T|, the molecular model for the fluid and the matrix are presented, and different 



aspects of the calculation are discussed. The results are given and discussed in section |V. 



II. THEORY 

In this section, we first summarize the theoretical expressions on which the present 
calculation is based. The reader is referred to the original papers for details ||ll[. We then 



present and discuss the extension of the formalism needed to treat the problem at hand. 

The replica-symmetric Ornstein-Zernike (RSOZ) equations relating the total pair correla- 
tion functions hapir) of a quenched- annealed mixture to the corresponding direct correlation 
functions Ca/sir) have been first obtained by Given and Stell and read: 

^00 = Coo + Po Coo ® ^00, (la) 

hio = cio + Po cio (8) hoo + pi Cc ® ho, (lb) 

hu = cii + Po cio ® hoi + pi Cc ® /ill + pi Cfe (g) he, (ic) 

hc = Cc + pi Cc ® he, (Id) 

where the indices and 1 refer to the matrix and fluid particles respectively, ® denotes 
a convolution, and the indices c and b denote the so-called connected and blocking (or 
disconnected) parts of hu and cn that obey 

hu = hb + he, (2) 

cii = Cb + Cc. (3) 



The replica OCT [|TT| is applicable to systems with pairwise additive interaction poten- 
tials and relies on the separation of these potentials Wa/3{r) into a (repulsive) part u^p{r) 
and a perturbative contribution Uapir). The properties of the so-called reference system, 
in which particles interact through the reference potentials, are assumed to be known, and 
the central quantities of the theory are then renormalized potentials denoted as Cajsir). 
Following Ref. [ITI, we quote here the equations obtained when only hard-sphere reference 
potentials with diameters aap are considered and when the matrix particles interact with 
each other only through such a hard-core interaction (i.e., Wooij) = u^Q^r)). 

The renormalized potentials can be defined through the replica optimized random-phase 
approximation (ORPA) closure which reads (/? = l/{kBT) where ks is Boltzmann's constant 
and T is the temperature) 

Coo('^) = c^o(^) for ^ > cToo, (4a) 

Cio(r) = cfo(r) - /?Mio(r) for r > aio, (4b) 

cii(r) = cfi(r) - /?Mii(r) for r > an, (4c) 

Cb{r) = cj^ir), (4d) 



complemented with the core conditions for the pair distribution functions 

fi'oo(^) = for r < o-qo, (5a) 

gio{r) = for r < a^, (5b) 

gii{r) = for r < an, (5c) 

(there is no core condition for gb{r) = hi,{r) + 1). Then, one has for the renormahzed 
potentials: 

Cio(r) = /iio(r)-/ifo(r), (6a) 

Cii(r) = /in(r)-/ifi(r), (6b) 

Cbir) = hir)-h^ir). (6c) 

The ORPA free energy of the quenched-annealed system can be determined in a closed 
form as 

— T H ^cx 

-AoRPA = 



V 
a"" + \[pl{cM - cf (0)) + 2pipo(aio(0) - 4(0))] - f [cn(0) - cf,(0)] 

rfk{ln[l - pi c,{k)] - Pl—-[c^{k) + p.-^^li'^l—]} 



(7) 



or, equivalently, by using the RSOZ equations, 

vAoRPA = 

a"" + ^[p?(ce(0) - cf (0)) + piPo(£io(0) - cfo(O))] - |[cn(0) - cf,(0)] 

c/klni^i^-^^/rfr,oi(r).oi(r) ^'^ 



2(27r3)y l-piC«(A;) 

+ ^JdrMr) - h^{r)]c^{r) + ^ J dr[ho,{r) - /.?i(r)]4(r), 

D 

where ^ is the free energy of the reference system. 

This ORPA free energy corresponds to the so-called energy route, since it satisfies the 
Gibbs-Helmholtz equation 
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(9) 
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where Uex is the configurational internal energy of the fluid inside the matrix and is given 
by 



-^ = ^ / dr[pl gnir) un{r) + 2pi po goiir) uoiir)]. (10) 



The mean-spherical approximation (MSA) is obtained from the ORPA when the reference 
system is described by the Percus-Yevick (FY) approximation, which corresponds to 

c^Q{r) = for r > aoo, (Ha) 

cfo(r) = for r > aio, (lib) 

cfi(r) = for r > (Til, (He) 

cf{r) = for r > 0. (lid) 

Within this latter approximation, an expression for the excess chemical potential can be 
derived: 

/3[/.ex - Z^^.] = ^[cii(O) - cfi(O)] - pi[Cii(0) - cfi(O)] 

-po[aoi(0)-c^i(0)] 

f3 f (12) 

= - / dr[pi gii{r) uii{r) + po goi{r) Uoi{r)] 

-y[cii(0)-cfi(0)]-^[coi(0)-c^i(0)]. 

As is often the case with simple closures, both the ORPA and the MSA are thermody- 
namically inconsistent, which means that the free energy obtained from the above energy 
route differs from that obtained by integration with respect to the fluid density of the com- 
pressibility equation |^,pO[ 

^ = l-piCe(0), (13) 

PiX 
where % is the isothermal compressibility of the fluid inside the matrix. This is an important 
remark, since it turns out that the compressibility route does not yield a critical point in 3 
dimensions [^, whereas the energy route provides one, albeit with an associated classical 
(mean-field) behavior. 

Finally, by truncating the OCT diagrammatic expansion to the first order beyond ORPA, 
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one obtains the 0RPA-I-B2 approximation for the free energy ||Tl[, that is the sum of the 
ORPA result, Eq. (9), and of the corrective term 

B2 = IpiPo I drhUr)[C,,{r)f + prpo j dv g^,{r)Y, -^[C,,{r)T 

^ ^ n=3 "" 

drhUr)[C^^{r)f + ^pl /rfr (?fi(r) f] l[Cii(r)]" (14) 

•^ n=3 ■ 

-\pI fdrh^{r)[C,{r)r-^pl f dr g^{r)Y,hc,{r)r. 

•^ -^ n=3 

The associated EXP approximation for the pair distribution functions is given by [^ 

^or''(r)=^o^i(r)exp[Coi(r)] (15a) 

^n''''(0 = ^n(r)exp[Cii(r)], (15b) 

g!''''{r) = g^{r)eMC,{r)l (15c) 

9r'{r) = g'^r{r)-g!^'{r). (15d) 



A noticeable property of the above equations describing the structure and the thermody- 
namics of a fluid adsorbed in a disordered matrix is that the quenched matrix enters solely 
through its pair correlation function /ioo(^)- The basic strategy followed in this paper is to 
take advantage of this property to study fluid adsorption in a model aerogel obtained by an 
out-of-equilibrium aggregation process (see below). The aerogel is then only characterized 
by its pair distribution function (in fact, its structure factor obtained by Fourier transform) 
determined by computer simulation. At first sight, this may lead to a conceptual difficulty. 
Indeed, the original treatment of quenched-annealed mixtures was restricted to situations in 
which the matrix is obtained by quenching an equilibrium distribution |^|-§]. However, as 
was stressed by Madden |^, there is no such difficulty and the RSOZ equations apply in 
a more general context, although the diagrammatic structure of the direct correlation func- 
tions is non-standard (indeed, these equations simply follow from the Legendre transform 
that relates the generating functional of the total correlation functions to the generating 
functional of the direct correlation functions). One can use the OCT approximation scheme 
as well. 

III. MODEL SYSTEMS AND NUMERICAL METHODS 

As stressed above, aerogels display very specific structural properties that can be ex- 
pected to have a strong influence on the properties of the adsorbed fluids. First, they are 
made of long chains of connected, almost spherical, colloidal Si02 particles. Despite the 
very low density of the adsorbent, this connectedness induces short-range correlations that 
differ from that generated by equilibrium hard-sphere matrices. Secondly, the silica strands 
forming the gel are organized into a long-range, fractally correlated structure whose experi- 
mental signature is a pronounced peak at small wave-vectors in the intensities measured by 
small-angle x-ray or neutron scattering. Accordingly, it seems important to use a model- 
ing of the aerogel structure that reproduces these two properties. In the present work, the 
aerogel is modeled as a connected cluster of spheres generated by a diffusion-limited cluster- 
cluster aggregation (DLCA) process |TB|,|TP|]. This process mimics the mechanisms involved 



in the synthesis of real aerogels, and the three-dimensional off-lattice version used here has 
been shown to provide gel structures that correctly account for both the short-range and 
long-range structural properties of real silica aerogels |1^ . 



For a given matrix density pq, the only information about the matrix that is needed in 
our calculation is the sample averaged structure factor Soo{q) defined from the averaged pair 
correlation function /loo(^) as 

5oo(g) = l+Po rhooirf-^^^^Anr'dr. 
Jo Qr 

Soo{q) has been computed following the lines of the original work of Hasmy et al. fl^ and the 



reader is referred to the corresponding publications for details. The corresponding curves 
for the three densities studied in this work (po = 0.025, 0.05, 0.1) are reported in Fig. |l|. 
Two characteristic features of aerogel structures are clearly seen on these data: at large q 
(gcToo > 27r), one finds for all densities 

sin(gcroo) 

600(g) - 1 + 2 , 

go-QO 



which corresponds to the fact that each aerogel particle has on average two neighbors con- 
nected to it, whereas at low q (go"oo < 2), a broad peak whose amplitude increases as density 
decreases is found, that originates from the long-range, fractal correlations of the particle 
strands forming the material. 

The matrix-fluid and fiuid-fiuid interactions are split into a hard-sphere part and a 
Lennard- Jones 12-6 tail following Weeks, Chandler and Andersen ||2^, i.e., Uais^r) = —Sas 
for a^p < r < l^^l^o^p and Uafiir) = 4ea^[(cra^/r)^2 - (aaf^/r)^] for r > 2'^/^aaf3- The use 
of hard-sphere cores for all the repulsive potentials is chosen to avoid difficulties associ- 
ated with calculating the free energy of fluids with soft-core repulsive forces. Moreover, the 
Lennard- Jones tails are truncated at r/aap = 2.5. 



In order to make contact with previous work |Tl|], we have only considered systems in 
which all particles have the same diameter, i.e., aoo = o"oi = (Tu = cr. It has to be stressed 
that this choice is rather unrealistic with respect to the actual experimental situation, since 
the silica particles forming the aerogels are always much larger than the adsorbate (typically, 
by a factor of ten). However, the model captures most of the specific nature of the correlated 
disorder induced by an aerogel and this is expected to be sufficient to give an account of the 
generic experimental behavior. Like in previous studies, we introduce the interaction ratio 
y = ^01 /^ii) and, in the following, a will be used as the length unit, eu as the energy unit 
and eu/ks as the temperature unit. 

In all calculations, the correlation functions of the reference system were obtained by 
solving the RSOZ equations in the Percus-Yevick approximation and the corresponding 
pressure and chemical potential were obtained by integration of the compressibility equation 
of state, Eq. (13). This is an additional approximation in the ORPA+B2/EXP framework. 
However, we know from previous studies |^ that the PY theory gives a good description of 
the reference-system properties when matrix and fluid particles have the same size. However, 
because these previous studies have dealt only with matrices quenched from an equilibrium 
distribution, we have tried to obtain another assessment of this approximation for the present 
situation of an "out-of-equilibrium" matrix by computing the Henry's constant Kh of the 
reference system. For a hard-sphere system, K^ is simply the fraction of the total volume 
of the matrix that is accessible to the center of a single adsorbate molecule. Thus, we 
have first estimated K^ by direct Monte-Carlo integration for various realizations of the 
model aerogels at the three working densities. Then, we have computed the theoretical 
value by using the (scaled-particle-theory) charging process expression (remembering that 
ir^ = exp[-/3/if,(pi = 0)]), 

Kh = exp ( -po / Goi{r; r; pi = 0)47rr^cir 

where Goi('^; r; pi = 0) is the contact value at zero adsorbate density of the fiuid-matrix pair 



distribution function when the corresponding hard-core diameter equals r p4[. The value 
of GQi{r;r; pi = 0) is determined in the Percus-Yevick approximation for various values 
of r between and a and the integral is estimated as a discrete sum. Note that a similar 



calculation has already been proposed by Ford et al. ^^ and appeared to be quite successful. 
The results of both estimations are reported in table | and show a very good agreement, 
thereby confirming the accuracy of the PY description of the reference system. 

In order to solve the RSOZ equations, we used Gillan's method [^ with a grid size 
of 0.02o" and 16384 points in most of the calculations. Once the correlation functions were 



obtained, we deduced the renormalized potentials and computed the free energy to the chosen 
approximation. In order to calculate the phase diagram of the fluid inside the matrix, series 
of adsorption isotherms /i vs pi were generated. When the curve had a loop, similar to those 
seen in the van der Waals theory of equilibrium fluids, the densities p'^^ and p'[ of the two 
coexisting phases were determined by solving the system of equations 

np.) = np[\ mp'i) = mpi), 

where P = — dAldV\ is the grand potential density (or thermodynamic pressure of the 

fluid i). 

In Ref . [|lT| , it was shown that there is a great sensitivity of the calculated phase diagrams 
to the approximations used in the theory. This is the case in the present calculation as well, 
as can be seen in Fig. 0, where the MSA and 0RPA+B2 phase diagrams are compared 
for typical values of the parameters. Note that the fluid density has been normalized by 
the porosity of the matrix that is equal to 1 — r^o, where r/o = rtp^a^ /Q is the corresponding 
packing fraction. In the following, we will concentrate on the results of the 0RPA+B2/EXP 
approximation. There are two motivations to this choice. The first one is of fundamental 
nature. Experiments show unambiguously that it is the low-density part of the phase dia- 
gram that is the most strongly affected by the presence of an aerogel. So, it is sensible to 
address the problem with a theory well adapted to this low-density regime. For bulk fluids, 
it is well known that, at variance with the lower-order approximations ORPA or MSA, the 
0RPA+B2 theory predicts the free energy with the same level of accuracy in the low-density 
gas and the high-density liquid phases and that the EXP approximation gives the exact low- 
density limit for the pair distribution functions 0]. This is not true anymore when one 



considers quenched- annealed system (for instance, the 0RPA+B2 free energy does not yield 
in general the exact second virial coefficient for the fluid inside the matrix, because the 
term of order p\ contains contributions from all orders in po), but one still expects that the 
0RPA+B2/EXP theory will remain accurate in the case of a low-density gas adsorbed in 
a low-density matrix. The second reason is more pragmatic and comes from the compar- 
ison with the simulation results of Alvarez et al. [^ for hard-sphere matrices of low and 



moderate densities. In this case, the 0RPA+B2 phase diagrams turn out to give the best 
agreement with the simulation results. This is shown in Fig. ^ where the simulation results 
tabulated in Ref. [^ are compared with the theoretical curves of Ref. 0] : the agreement 
is reasonably good (much better indeed that with MSA). 



IV. RESULTS 

We first present the results for the structure of the adsorbed fluid, since it shows remark- 
able features originating in the specific properties of the aerogel, in particular its connect- 
edness. 

The first effect of the aerogel structure shows up already in the case of purely repulsive 
matrix-fluid interactions {y = 0). This is illustrated in Fig. ^, where the matrix- fluid 
correlation function 5'oi(^) at gas density (pi = 0.02) displays a very clear depletion near 
contact. This "correlation hole" comes from the presence of a region surrounding each 
matrix particle from which the adsorbate molecules are excluded. This region disappears 
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when density is increased, as can be seen from the goi{r) at hquid density (pi = 0.70). This 
correlation hole, as well as the accompanying cusp at r/a = 2, is a well-known feature of 
molecular systems made of tangentially jointed hard-spheres, like our model aerogel. For 
instance, it has been found in fluids of freely-jointed hard-sphere chains, by both computer 
simulation and theory |2^ . It is a simple steric effect due to the presence of the neighboring 
beads that hinders the approach of a gas molecule. 

The second effect, expected to manifest itself at low enough temperatures when attractive 
matrix-fluid interactions are turned on, is a much stronger binding of the adsorbate particles 
to the solid matrix than in the previously investigated models. Indeed, on one hand, two 
tangent aerogel beads form adsorption sites with potential energy equal or lower than — 2£:oi 
(the direct contribution of the two beads being still lowered by contributions from their 
neighboring beads); on the other hand, the aerogel structure presents large voids [0,0 in 
which the adsorption potential energy due to the matrix is zero. This leads to a highly inho- 
mogeneous potential energy field and thus to strong preferential adsorption on the aerogel 
strands. The situation is very different in the case of an equilibrium hard-sphere matrix 
at a similar density, as shown by the following simple calculation. For such a matrix at a 
density po = 0.05, the average distance between two spheres is roughly d = {po)~^^^ — 2.7. 
If an adsorbate particle is equidistant from two matrix particles, it experiences from these 
particles a potential energy roughly equal to 2uio{r = d/2) = —1.075 Eqi, whereas if it is 
adsorbed on a matrix particle, it experiences a potential energy simply equal to — ^oi- The 
situation is thus far less contrasted than in an aerogel and one expects a relatively weak 
binding of the adsorbate on the matrix. 

This second effect should result in highly inhomogeneous distributions of the fluid parti- 
cles and is thus difficult to detect from averaged quantities like pair distribution functions. 
We find nevertheless a few indications of its infiuence on the adsorbate structure. One of 
them can be seen in Fig. ^. It is indeed quite direct to see from the values of the pair dis- 
tribution functions near contact that, at the two vapor densities considered (pi = 0.02 and 
pi = 0.2), there are roughly twice more adsorbate particles around a matrix particle than 
around a given adsorbate particle. This is compatible with the picture of the formation of 
a strongly bound monolayer of particles adsorbed on the aerogel strands. It agrees with the 
statement by Chan et al. that the "gas" phase in very dilute aerogels actually consists of va- 
por plus a liquid film that coats the random silica strands p. This situation differs strongly 
from that in equilibrium hard-sphere matrices, where the prevailing adsorption mechanism 
is expected to be the formation of liquids droplets in the denser regions of the matrix. In 
this latter case, the environments of the matrix and fluid particles are not too different. This 
is precisely what is found in Ref. [O (see Fig. 13 of this paper and the related discussion). 



Another indication of strong adsorption on the aerogel strands is the presence of structural 
correlations typical of the aerogel structure in the adsorbed low-density phase. Indeed, if the 
adsorbate particles coat the aerogel strands, it is expected that their distribution will retain 
part of the spatial arrangement of the substrate, in particular the long-range correlations 
that show up in the form of a maximum at small wave vectors. To test this idea, we have 



calculated the partial structure factors of the system defined as [^ 



Po + Pi 

Po + Pi 

Results are shown in Fig. H, where it is clearly seen that a broad correlation peak builds up 
at small q {qa ~ 0.2) when increasing the vapor density, at a position in good agreement 
with the position of the corresponding peak in the aerogel structure factor (compare with 
Fig. 1). This can again be interpreted as evidence for the progressive formation of a thin 
layer of adsorbed particles coating the gel strands. 

These characteristic features of adsorbed fluid phases in aerogels are important at low 
adsorbate density only. Indeed, as can be seen in Figs. || to ||, they are almost completely 
absent at liquid densities, for which the physics of the system is dominated by the hard- 
core interaction between adsorbate particles so that the structure of the adsorbent becomes 
of marginal importance. Thus, these observations stress a posteriori the need to use a 
theoretical framework well adapted to the description of the low-density regime, hence the 
use of the 0RPA-I-B2/EXP theory in this work. But they also show that the case under 
study is an especially difficult one in the context of the OCT: indeed, of the two effects 
discussed above, the first one, that is purely steric in origin, competes with the second one, 
that is only due to the attractive matrix-fiuid interactions. 

We now turn to the presentation of the capillary phase diagrams. We only report the re- 
sults corresponding to the case where both fiuid-fiuid and matrix-fiuid attractive interactions 
are present, which is the case of experimental relevance. 

Fig.0 shows coexistence envelopes of the fiuid phases obtained from the 0RPA+B2 ap- 
proximation as one varies the ratio y for a typical aerogel density po = 0.05 corresponding 
to a porosity 1 — tjq = 97.4%. It is seen that the variation of the capillary critical tem- 
perature with y is non-monotonic: Tec first increases and then decreases as one increases 
y, the maximum value being found around y = 1. On the other hand, the critical density 
increases monotonically and becomes larger than the bulk value for y = I (this not due to 
the normalization by 1 — rjo). The latter behavior agrees with the experimental observation 
that pic is displaced towards the liquid phase when the interaction with the porous material 
is sufficiently attractive 0. Despite the low density of the aerogel, there is also a significant 
narrowing of the coexistence curves compared to the bulk one when y > 1, although this 
effect is still smaller than that observed in the experiments (see below) 0,^]. 

Fig. 1^ shows the effect of varying the aerogel density at a fixed value of the interaction 
ratio y = 1.5. It is found that the normalized critical density increases only slightly with 
the matrix density and that the width of the coexistence curve does not vary very much. 
However, for y = 1.5 and po = 0.025, a "precondensation" transition is obtained, similar 
to what has already been found in various models @,0] (see also Maritan et al. |^1[). A 
surprising result seen in this figure is the non-monotonic variation of the capillary critical 
temperature that first increases with po and then decreases, leading to the somewhat un- 
expected result that for a low enough matrix density Tec is larger than its bulk value. The 
physical origin, if any, of this behavior is unclear to us, but there is a priori no fundamental 
reason forbidding such a behavior. A possible explanation could be the very weak confine- 
ment effect due to the dilute gel. Indeed, it is well known, already at the mean-field level, 
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that confinement induces a lowering of tlie critical temperature of a fluid. In particular, this 
is what is found in all cases of simple pore geometries [Q . But aerogels are rather special in 
this respect: because of their very open structure, they have no real pores with well defined 
walls. The comparison with the case of equilibrium hard-sphere matrices gives an insight 
on the potential relevance of this fact. Indeed, in the latter case, it has been found that 
adsorption partly counterbalances the influence of confinement, so that the lowering of the 
critical temperatures is less pronounced in the case of attractive fluid-matrix interactions 
than in the completely repulsive case. A possible origin of our finding for aerogels could thus 
be that the effect of the attractive matrix- fluid interaction now dominates the very weak 
influence of confinement. 

The remaining question is of course to see how our results based on a realistic description 
of the structural correlations in an aerogel compare with experiments and with previous 
results obtained with the less realistic equilibrium hard-sphere model. 

As already stated in the introduction, the most significant experimental observation is 
a marked narrowing of the liquid-vapor coexistence curve relative to that in the bulk, a 
narrowing that is accompanied by a slight displacement of the critical point to a lower 
temperature and a higher density. Chan et al. [0 have proposed to quantify this narrowing 
by fitting the top of the coexistence curves for ^He and N2 adsorbed in 95% aerogels with 
the critical expression 



Pi- Pv ^ Q 
Pic 



Tec - T 
T 



where p^ and pi are, respectively, the vapor and liquid densities of the adsorbed fluid, and 
they have found that the amplitude factors B are respectively about 14 and 2.6 times 
smaller than the bulk values for these two fluids (in both cases, the value of the exponent 
j3 is consistent with that found for the bulk fluids). These values, which seem indicative 
of a large reduction of the phase coexistence envelopes, have nevertheless to be taken with 
care. Indeed, in the case of '^He, the previous result has been obtained from data in the 
temperature range 5.15-5.17 K, T^c being found equal to 5.167 K. Further experiments in a 
wider temperature domain (2.5-5 K) have shown that the reduction of the phase envelope 
is much more pronounced near the critical point than well below the critical temperature 
||32|| . For instance, when the temperature is roughly 0.64 Tc, the ratio of the widths of the 
bulk and adsorbed fluid phase diagrams is not larger than about 2.2. 

As can be seen in Figs. ^ and |, for values of y that are large enough, one finds that our 
model gives results in good qualitative agreement with the picture emerging from experi- 
ments: the critical density is slightly shifted to values higher than in the bulk and, except for 
the lowest aerogel density, the critical temperature is found to decrease, albeit moderately. 
A significant narrowing of the coexistence curves is also found. For instance, with po = 0.05 
and y = 1.5, the ratio of the widths of the bulk and adsorbed fluid phase diagrams is about 
1.4 at T = 0.8 ~ 0.65 Tc and it increases with y. (However, our present approach is unable 
to capture the behavior seen in the "^He experiment in the close vicinity of Tec-) 

Figs. H (a) and ^ (b) show a comparison between the coexistence curves obtained for a 
fluid adsorbed in model aerogels and equilibrium hard-sphere matrices. It is seen that, at 
given matrix density po, and ratio of interactions y, the coexistence curves in the case of the 
model aerogel are always narrower, and the critical temperatures always higher, than in the 
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equilibrium hard-sphere case. The effect is quite significant for y = 1.5, where the ratio of 
the widths of the bulk and adsorbed fluid phase envelopes is about 90 % at T = 0.8 ~ 0.65 T^ 
for hard-sphere matrices, compared to about 70 % for aerogels. Combined with the fact that 
the width of the coexistence curves varies only weakly with the aerogel density, we thus find 
that using a more realistic model of the adsorbent makes it possible to obtain a behavior 
that more closely resembles the experimental one, i.e., a substantial narrowing of the phase 
envelope associated to weak changes in the critical temperature and density. (It should 
also be stressed that for the small matrix densities used here, the hard-sphere matrix is far 
from being a connected medium.) We note however that this improvement of the situation 
compared to experiments is limited to a certain range of parameters, since the differences 
between the two types of matrices become rather modest for larger values oi y {y > 2), for 
which it seems that the strength of the attractive matrix-fluid potential limits the relevance 
of the detailed structure of the adsorbent. 



V. CONCLUSION 

In this paper, we have presented a theoretical study of the phase diagram and the struc- 
ture of simple fluids adsorbed in high-porosity aerogels. For these systems, the perturbation 
induced by the adsorbent is a priori small, and one expects the integral-equation approach 
combined with the replica formalism to provide a reasonable description. The novelty of the 
present work is that a realistic description of the aerogel structure has been incorporated 
via the use of the static structure factor of model aerogels built from an off-lattice diffusion 
limited cluster-cluster aggregation process. Unlike cruder models based on quenched config- 
urations of equilibrium hard or penetrable spheres, the correlations between matrix particles 
at the low density studied (corresponding to porosities between 90% and 97.5%) keep track 
of the short-range connectedness of the aerogel and of the long-range fractal behavior. The 
predictions of the theory are in qualitative agreement with the experimental results, show- 
ing a significant narrowing of the gas-liquid coexistence curve associated with weak changes 
in the critical temperature and density. In addition, the influence of the aerogel structure 
is shown to be important at low fluid densities. Improving the present approach is still 
needed in two directions: first, to provide a better description of the critical region, but 
this represents a very challenging task for systems like fluids adsorbed in disordered porous 
media that are akin to random fields models; secondly, to go one step further toward a 
realistic modeling of the adsorbent by allowing the beads composing the silica strands of the 
aerogel to be much larger (typically, by a factor of ten) than the fluid particles. The latter 
improvement of the description requires both the development of more appropriate closures 
of the Ornstein-Zernike equations for mixtures of particles with very different sizes and the 
introduction of a more realistic potential between the matrix and the fluid particles (using 
for instance the composite-sphere potential proposed in Ref. [p3|). 
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TABLES 



po MC PY EHS 

0.025 0.934 0.927 0.898 

0.05 0.868 0.856 0.802 

0.1 0.739 0.718 0.629 



TABLE I. Henry's constant of a hard-sphere fluid adsorbed m DLCA aerogels, calculated by 
du'ect Monte-Carlo integration (MC) and within the Percus-Yevick approximation (PY) through 
a scaled-particle-theory charging process. The fluid and matrix particles have the same diameter. 
For comparison, the values for hard-sphere matrices (EHS) calculated as in Ref. |11| are also given. 
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FIGURES 
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FIG. 1. DLCA aerogel structure factors at three different densities: po = 0.025, 0.05, 0.1 (from 
top to bottom). These curves represent averages over 40 realizations of the aerogel generated at a 
given density. 
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FIG. 2. Comparison of the coexistence envelopes for a fluid confined in an aerogel matrix 
(po = 0.05, y = 1.5) computed within the MSA (solid line) and the 0RPA+B2 (dashed line). 




FIG. 3. Comparison of the coexistence envelopes for a fluid confined in purely repulsive (y = 0) 
equilibrium hard-sphere matrices (po = (bulk), 0.046 and 0.15 from top to bottom) obtained by 
computer simulations (symbols) p^ and within the 0RPA-I-B2 (lines) pi]. 
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FIG. 4. Fluid-fluid distribution function gii{r) and matrix-fluid distribution function goi{r) 

calculated with the EXP approximation in an aerogel matrix with po = 0.05 and y = at T = 0.9. 

Solid line: pi = 0.02; dashed line: pi = 0.70. 
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FIG. 5. Fluid-fluid distribution function gii{r) and matrix-fluid distribution function goi{r) 
calculated with the EXP approximation in an aerogel matrix with po = 0.05 and y = 1.5 at T = 0.9. 
Solid line: pi = 0.02; dashed line: pi = 0.2; dotted line: pi = 0.72. 
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FIG. 6. Partial fluid-fluid structure factor Sfj^{q) and matrix-fluid structure factor SQ^{q) 
calculated within the EXP approximation in an aerogel matrix with po = 0.05 and y = 1.5 at 
T = 0.9. Solid line: pi = 0.02; dashed line: pi = 0.2; dotted line: pi = 0.72. Note that the 
oscillations of Sq-^^^q) around zero at large q prevent us from using a logarithmic scale for Sq-^^Iq). 
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FIG. 7. 0RPA+B2 predictions for the coexistence envelopes of the fluid phases in an aerogel 
matrix with density po = 0.05. From left to right: y = 0.5, 1.0, 1.5, 2.0, 2.5. The solid line is for 
the bulk fluid. 



23 




0.4 0.6 

Pl/{1-Tl0) 



FIG. 8. 0RPA+B2 predictions for the coexistence envelopes of the fluid phases in an aerogel 
matrix with interaction ratio y = 1.5. From top to bottom, po = 0.025, 0.05, 0.1. The solid line is 
for the bulk fluid. 
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FIG. 9. Comparison of the gas-liquid coexistence envelopes predicted by 0RPA-I-B2 in an 
equilibrium hard-sphere matrix and a DLCA aerogel at the same matrix density po = 0.05. (a) 
y = 1.5. Dotted line: hard-sphere matrix; dashed line: aerogel, (b) y = 2. Dotted line: hard-sphere 
matrix; dashed line: aerogel. The solid line is for the bulk fluid. 
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